Bose and Fermi gases in the early universe with self-gravitational effect 
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We study the self-gravitational efTect on the equation of state (EoS) of Bose and Fermi gases in 
thermal equilibrium at the end of reheating, the period after quark-hadron transition and before 
Big Bang Nucleosynthesis (BBN). After introducing new grand canonical partition functions based 
on the work of Uhlenbeck and Gropper, we notice some interesting features of the newly developed 
EoSs with distinct behaviors of relativistic and non-relativistic gases under self-gravity. The usual 
negligence of the self-gravitational effect when solving the background expansion of the early universe 
is justified with numerical results, showing the magnitude of the self-gravitational modification of the 
state constant to be less than O{10~^^). This helps us to clarify the background thermal evolution 
of the primordial patch. Such clarification is crucial in testing gravity theories, evaluating infiation 
models and determining element abundances in BBN. 

PACS numbers: 98.80.-k, 04.40.-b, 05.30.-d 



I. INTRODUCTION 



In most history of our universe, gravity has been the 
dominant force in determining cosmological evolutions. 
Due to the onset of gravitational instability, primordial 
fluctuations generated in cosmic inflation start growing 
by attracting dark matter and gases nearby after the 
matter radiation equilibrium at a temperature approx- 
imately 4000 K m . Consequently, structures of the first 
generation galaxies are seeded. Shortly after the forming 
of proto-galaxies, the hydrogen and helium gases within 
them begin to condensate and make the first stars. In 
the formation of any astrophysical body, self-gravity is 
beyond question one of the key factors in determining 
the specific destinies of stars, quasars, black holes, and 
neutron stars. For example, the self-gravity of neutron 
stars is considered using the affine model 2, 3] , where the 
effective relativistic self-gravitational potential is approx- 
imated with Tolman-Oppenheimer-Volkofl[ (TOV) stel- 
lar structure equations. Indeed, gravity, though much 
weaker than the other three fundamental forces, i.e., 
strong, weak and electromagnetic forces, with its long 
range of action and omnipresent character from equiva- 
lence principle of general relativity, determines most as- 
tronomical evolutions. One persuasive example showing 
the significance of self-gravity is the boson star. Though 
massive boson stars are prevented from the Heisenberg 
uncertainty principle, they can actually be formed in the 
early universe when the self-interaction of scalar par- 
ticles is considered d, Various methods are sub- 
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sequently developed to study different features of bo- 
son stars. People deal with a system of self-gravitating 
bosons either by solving the Schrodinger and Poisson 
equations in a Newtonian way 0; or by fully account- 
ing the relativistic effect using Klein-Gordon and Ein- 
stein equations Jj. When the Einstein gravity and the 
Klein- Gordon field couple with each other, exotic bosons 
which possess a small mass may undergo Bose-Einstein 
condensation driven by gravity and collapse to form bo- 
son stars in the early universe or latter in galaxy center 
as a candidate for dark matter (sl-fioj. 

If self-gravity can trigger phenomena at the early uni- 
verse so dramatic like boson stars or primordial black 
holes [ill, [13 ^ot mentioned before, the self-gravitational 
effect in the primordial patch of the early universe should 
not be easily discarded without a robust proof. An exact 
evaluation of the self-gravity of relativistic particles at 
the extremely high energy scale before Big Bang Nucle- 
osynthesis (BBN) (not high enough to reach the boson 
star critical temperature) becomes unavoidable. How- 
ever, most discussions of boson stars are at relatively 
high energy scale ~ 1 GeV [13], and at the same time 
entail the non-zero mass of scalar particles to be greater 
than around 10 eV to preserve radial stability [lj|. But 
during the radiation dominated era before BBN, main 
ingredients are massless photons, and the homogene- 
ity of the early patch is preserved unlike that in boson 
stars. As a result, it is difficult to fully perceive the self- 
gravitational behavior of photons in the early universe 
by simply extending the methods developed for boson 
stars. Although there are several attempts to include the 
self-gravitational effect of non-relativistic gases at equi- 
librium, such as Monte Carlo simulations, analytic mean 
field methods and low density expansions [l5| , no method 



2 



so far exists to directly solve the self-gravity of massless 
particles, i.e., photons. As a result, we try to develop 
an appropriate way to deal with the self-gravity of mass- 
less bosons in the early universe. The newly developed 
method will be applied to massive particles as well. 

Before BBN, the universe can be well described as in 
thermal equilibrium with very high temperature and fre- 
quent interactions. Remarkable features of the early uni- 
verse, homogeneity, isotropy and thermal equilibrium to- 
gether can help us simplify the discussion to a great ex- 
tent. Since the Equation of State (EoS) can serve as a 
good description of a homogeneous systems, we manage 
to delve into the problem with statistical physics. As 
is well known, once including self-interactions between 
particles, the corresponding partition function shall be 
changed. Such change further leads to modifications in 
the EoS. We try to uncover the self-gravitational effect 
by examining the corresponding statistical properties of 
Bose and Fermi gases first. Whether corrections from 
self-gravity are negligible or bear some new features not 
considered before shall be tested in detail in this work. 

In the radiation-dominated epoch, gravities between 
relativistic particles are commonly neglected in existing 
models. When the temperature of the universe decreases 
to approximately 10^^ K, particles are in thermal equi- 
librium without strong nuclear interactions and the prin- 
cipal ingredients at that time are photons, neutrinos, 
electrons and their anti-particles. With the energy of 
the universe dominated by photons, other ingredients, 
whether massive particles like electrons and positrons or 
the massless ones, all behave radiation-like due to their 
frequent collisions with photons and each other at the 
ultra-high temperature. The matter soup at that time is 
suggested to follow the equation of state P = yo/3 with- 
out considering gravitational attractions between parti- 
cles. From the equivalence principle, massless particles 
with both high kinetic energies and particle densities 



where a — —fi/ksT represents the chemical potential, 
and m is the mass of the particle. 

A. The non-relativistic case 

For non-relativistic particles, which satisfy the disper- 
sion relation Ek = p^/2m, the first integral in Eq. ([2]) 



may possess a considerable gravitational effect. We deal 
with self-gravitational interactions between these ther- 
malized particles in a Newtonian way. With a statisti- 
cal approach, we can avoid the theoretical difficulties in 
solving equations of general relativity as well as numer- 
ical errors in actual realization, and at the same time, 
perceive macro-behaviors without any impediment from 
micro-complexity. In Sec. |IT1 we deduce a new form of 
grand canonical function for Bose and Fermi gases in the 
relativistic and the non-relativistic limits. In Sec. IIIH we 
apply the newly-developed partition function to the early 
universe before BBN. We also analyze the time depen- 
dent term of the EoS and set a strict upper bound on the 
modifications of the state constant from self-gravity. Pos- 
sible influences on relating topics are discussed in Sec. lIVI 

II. GRAND CANONICAL PARTITION 
FUNCTION FOR SELF-GRAVITATIONAL BOSE 
AND FERMI GASES 

In standard statistical physics, people deal with inter- 
actions between gases by including e~'^^^/^'^ in the grand 
canonical partition function, where (pij is the potential 
energy resulting from the interactions between two par- 
ticles i, j. The potential energies (pij are set to be zero 
in general. Such approximation works pretty well when 
interaction is negligible, that is, the density-temperature 
ratio is extremely low. On the other hand, if the density 
is relatively high like that in the early universe, the factor 
from the self-gravitational effect of bosons and fermions 
becomes appreciable. As Uhlenbeck and Cropper [l^ 
proved, the partition function of non-ideal Bose or Fermi 
gases should be modified as 

g-0„/feT _^ [-AnmkBTrlj/h^]) , (1) 

where the plus sign is for Fermi gas. 

The grand canonical partition function becomes 



(2) 



I 

turns 



/ e ^^^^^dpi = {2mT:kBT)~ . (3) 

i 

Since in the non-relativistic case, 0.^ = Grriimj /rfj, the 
second integral becomes 



-aN r E. r'^<3 



^ - E mhm n J ' dp, Jl[e-^BT(i± exp [-AnmkBTrl/h'] ) dr. 
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(1 ± exp [-ATTmkBTrfj/h^])dr, 



1 ± ^ exp [-AnmkBTrl/h^] dr^ . (4) 



r 



With exp [-47rmfcsTry/i2] < 1, the omissions made 
in Eq. ^ are justified. At the same time, we use mean 
field theory to include only the average effect of r,j . The 
thermal equilibrium used in mean field theory is guaran- 
teed in the early universe. As — > oo we approximate 
N{N - l)/2 ^ following which Eq. g]) becomes 

e i-r.kgT \ I ± — [-47rTOfcsTrJ//i2] j dr*dr 
V^'^ I I e'^^'^^^dr, ± / —exp[-ATTmkBTrl/h'^] 



This reveals that the upper bound of the integration cor- 
responds to the spacial scale of the system. 



The second term of Eq. ^ is 
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e ^r.kgT _ 1 exp [-AnmkBTrl/h^] dr 



where is the average distance between two parti- 
cles. We can see that the first term on the right side 
of Eq. ([5]) represents the gravitational correction of the 
grand canonical partition function in conventional treat- 
ment and the second term corresponds to usual Bose- 
Einstein or Fermi-Dirac statistics, while the last term 
is a newly-introduced coupling term from the gravity of 
bosons or fermions. Again, the upper sign is for Fermi 
gas. 

In the early universe, particle number and energy den- 
sity are relatively high with temperature^ 10^^ K, i.e., 
exTp[—N^Gm'^ /2r^,kBT] ^ 1. To integrate the first term 
of Eq. (O, we expand it into the the first-order Taylor 
series and get 
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2r^,kBT ^ 
16 / kBT 



(6) 



dr* 

(5) 



j —exp [-^TimkBTrl/h'^] dr, 

4:71 —exp[~4:7rmkBTrl/h^]rldr^ 
Jo 2 

f°° N"^ 

4tt / —exp [-4nmkBTrl/h'^] r^dr. 
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(7) 



where Anonrcia = h/y^2'KmkBT is the thermal de Broglie 
wavelength, which should be much smaller than the ra- 
dius of our interested system. If the scale of the system 
is equal to or smaller than the thermal de Broglie wave- 
length, a situation which might happen at a much ear- 
lier time, such approximation may not be valid, and we 
should use the error function erf (x) in the integration in- 
stead. Again applying the first-order Taylor expansion, 
we get the remaining term 
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/ exp [—AnmkBTr'^ / h^] r,dr, = 



N^TrGm^ 
4kBT 



nonrcla' 



(8) 



Finally, we find the new grand canonical partition func- 
tion 
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nonrcla 



(9) 
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with which we can obtain the total energy hy E = 
dlnE/d(3 with /3 = l/fc^T, and then the entropy. As sug- 
gested above, this partition function is apphcable when 
the scale of the system is much larger than the ther- 
mal de Broglie wavelength. However, such premises may 
be wavered at a much earlier stage. The first term of 
the partition function is the grand canonical function for 
ideal Boltzmann gas. The second term represents the 
self-gravitational correction of Boltzmann gas, and its 
magnitude corresponds to the portion that gravitational 
energy Ep ^ Gm? /V'^/^ contributes to the total energy 
^'totai ^ ksT. In turn it manifests whether gravitational 
correction is important in dealing with the evolution of 
the universe. The third term is the correction from iden- 
tical principle, and its magnitude (proportional to the 



cubic square of thermal de Broglie wavelength over the 
volume of the system V) evaluates how strong quantum 
identity is. As a whole, the self-gravitational effect on 
fermions is restricted by both the strength of gravity and 
the intensity of quantum effect. 



B. The relativistic case 

For highly relativistic particles, the rest mass is negli- 
gible compared with the kinetic energy, thus the disper- 
sion relation becomes = pc. Then the grand canonical 
partition function is found to be 



- ^ / e"^' ^dp, / e r^^'-BTc- 1 ± ^exp [~A-K'^p,kBTrl/h'^] dr,^-. (10) 



AT " •' \ i<j 



We can solve it exactly by diagonalizing X]i<j r- kllrc^ 
first. But such treatment is proved to be unnecessary 
since, on one hand, the self-gravitational effect in the 
relativistic case acts as a correction term, and on the 
other hand, for our problem in the homogeneous and 
isotropic Robertson- Walker metric background at excep- 
tionally high temperature, relativistic particles are in a 
thermal equilibrium state evolving in a radiation man- 
ner resulting from frequent collisions. Therefore, we ex- 
pect self-gravity alone will not significantly destroy this 
homogeneity directly shown by CMB observations. All 



particles in the patch then possess the same characteris- 
tic average momentum due to the frequent collision. We 
then approximate the summation as 

where p is the average momentum of the particles in the 
system. So the grand canonical partition function be- 
comes 



N ' i •' 



±n / e^V^+^^^^^+^^-'^^^r-dp^drJ 



/ (8^)"^ dr 

j^Nlh^^J GP Y 

\kBT ^ 2r,,c^kBT ) 

±E]^/7 ^ (11) 

AT J ^ l-V 47r2^^-— 



To integrate the first term in Eq. (jlip . we perform the tude of gravity being small outside the critical sphere 
Taylor expansion to the first order knowing the magni- with radius 
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\kBT ^ l^j 2ri,c^kBT 



3N 



{Stt) 



N 



s SAT 



dr,; 



(12) 



The integration from to r,„ — Gp/ (? 



Gp 



3N 



dr,; 



(13) 



is neglected since its magnitude is proportional to r"^^^, 
much smaller than the first one in our analysis when grav- 
ity contributes as a correction term. 

Under the same mean field approximation, Eq. (fT2)) 
becomes 



{ksr) 



3N \ ^ V 6 



(14) 



Before going on to integrate the second term of 
Eq. ([Tl]) . we first compare the magnitude of differ- 
ent terms in the denominator + J2j 2r 'c'^ksT ~^ 



,kBTrf 



„3/2 



E,4^'-^- When ^ ^4^2 ^^^^ , _ 

the corresponding energy density can be approximated as 
U ^ ^ 10*6 in the Planck unit. However, the actual 

energy density is of magnitude 10^^ [T^, meaning that 

Be- 



ksTrj 



the first term is 0(10^ ) compared with tt 
sides, since in our interested time period, the temperature 
varies within one order in magnitude, so long as the num- 
ber density does not change beyond 0(10^^), i.e. the vol- 
ume of the universe increases by less than 0(10^^) within 
Is (such circumstances are fully guaranteed in a classical 
radiation-dominated universe), -^^^ can then be safely 
neglected. But for some earlier stages this treatment 
may fail. In the same way, ^ .^ksT 0(10"''^) com- 
pared with the third term and can be justifiably omitted 
if the number density varies within 0(10'*°). Neglecting 
the first two terms after careful magnitude evaluation, 
we have 
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(8^) 
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3N 
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(15) 



Therefore only the first integration contributes to our grand canonical partition function for relativistic parti- 
final result, corresponding to a unification of Bose and cles including the self-gravitational effect becomes 
Fermi distributions under the relativistic limit. The 
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(16) 



where Aj-cia = ch/2i:^l^kBT is the thermal de Broglie 
wavelength for massless particles. In Eq. ([T5)) . the first 
term is merely the ordinary ideal Boltzmann gas term, 
while the second is induced by self-gravity, proportional 
to Gh/ K^ai&c^V^/^ , i.e. the gravitational potential energy 
over the self-energy of particles. By comparing the two 
terms we can catch a glimpse of how gravity acts on the 
statistical properties of the relativistic particles. 



III. RESULT IN THE EARLY UNIVERSE 



Before applying the newly developed method to a spe- 
cific background, we first review the whole evolution of 
the early universe. In the very early stages of the evolu- 
tion, to resolve several problems in the Big Bang cosmol- 
ogy, like the horizon problem, magnetic-monopole prob- 
lem, structure formation problem and so on, the uni- 
verse is suggested to undergo an infiation between time 
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scales 10~ s and 10~ s. By the end of the infla- 
tion, the background temperature is roughly 10^* K when 
strong, weak and electromagnetic forces are unified. As 
the temperature cools down, the strong force freezes out, 
weak and electromagnetic interactions are subsequently 
off stage and before 10""* s. Afterward, neutrinos decou- 
ple at T 10^ K. Then Big Bang Nucleosynthesis be- 
gins at temperature 10^ K 10^ K. At the end of BBN, 
photons decouple from matter, and form the cosmic mi- 
crowave background radiation we observe today. To use 
our statistical method in studying the self-gravitational 
effect, assumptions of homogeneity, isotropy and thermal 
equilibrium should be justified. We thus choose to study 
the universe before BBN when three criteria are largely 
met and the background temperature Tth is of magnitude 
10^^ K [17|, well below the corresponding thermal en- 



Their grand canonical partition function is 



ergy of the rest mass of tt mesons and strong interactions 
can be neglected. Matter and radiation in the universe 
at that time are in thermal equilibrium, resulting from 
the rapid collisions in the exceedingly high temperature 
despite of the fast universal expansion. By identifying 
particles whose rest masses are below the correspond- 
ing thermal energy of the background temperature Tth, 
we find the main ingredients to be electrons, neutrinos, 
their anti-particles and photons. These abundant genres 
are highly relativistic. For example, the rest energy of 
the electron corresponds to the thermal energy at tem- 
perature ^ 10^ K, two orders of magnitude smaller than 
the background temperature. And non-relativistic par- 
ticles, mainly protons and neutrons, constitute only an 
extremely small portion of the whole ingredients, approx- 
imately lO"""^" as many as that of relativistic ones. Be- 
low we deal with relativistic and non-relativistic particles 
separately, though the latter is unable to significantly in- 
fiuence the background evolution. 



A. Thermalized relativistic particles 



At a temperature higher than 8000 K, the universe is 
radiation dominated. Consequently, relativistic particles 
are undoubtedly the main ingredients during our study- 
ing period with a temperature approximately 10^^ K. 
And both massive and massless components act in a ra- 
diationlike manner due to frequent collisions. We use the 
Planck unit fcB = /i = c = G' = linall calculations 
and some of our discussions below. Moreover, the vol- 
ume Vinit , which is about four light years in radius in our 
partition function, corresponds to the biggest thermal- 
ized region constrained by CMB anisotropy and grows as 
the extension of particle horizon of the primordial patch. 
Though the exact size of the early universe at this time 
is not determined yet. The time span of our study is 
0.1 s ~ 1 s. We consider the relativistic bosons first 
and later find that fermions behave similarly as expected. 
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(17) 

The chemical potential is set to be zero since the particle 
number of the photon, which is the main ingredient, does 
not conserve. Here, instead of directly solving the Fried- 
man equation with the density and pressure derived in 
the new grand canonical partition function, we approxi- 
mate the temperature evolution the same as that in the 
classical radiation-dominated model neglecting the self- 
gravitational effect 



(18) 



where a is a constant to be determined by the initial 
condition. Combining Eqs. (IT71) and ^TE\\ . we numerically 
find that the density of relativistic particles relies on time 
as 



P - 



bi , h 



Vtv' 



(19) 



where bi and 61 are separately 6.58 * 10^^*^ J • s^/m^ and 
7.09 * 10^® J • s^. This differs from the conventional result 
without considering the self-gravitational effect where 
p ^ ^ That is, once including the self-gravity, 

energy density decreases much slower. This disparity in 
energy density can be explained as an increase in gravita- 
tional potential caused by the rapid elongation of particle 
separation in cosmic expansion. This remedies the con- 
ventional energy lost with the decrease of particle density 
and the extension of particle wavelength. The pressure 
behaves as 



c2 



(20) 



where 62 = 1-67 * 10 J • s . As a result, the state con- 
stant no longer stays invariant but instead relies on time 
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(21) 



where di has a numerical value 3.97 x 10^" in Planck 
unit which is 9.69 x 10~^^ s^/^ .ni^. Knowing that during 
our studying period t ^ lO'*^ jn Planck unit, or 0.1s in 
standard unit, and the volume is V ^ 10^^^ in Planck 
unit, or corresponds to four light years in radius justified 
by CMB observation, we find that the correction term 
in Eq. ^ is of the magnitude 0(10""'^). Since the 
universe is radiation dominated at this time, we find the 
correction term di/\/tV ^ t~^. Differing from the non- 
relativistic case below, the modification from self-gravity 
rapidly diminishes as time goes on and volume expands. 

Besides, we can roughly estimate the modification term 
at much earlier times, such as the end of inflation when 
t - 10-34 s and y ~ to be 0(10"^^). Thus 

the estimated magnitude of the modification is compara- 
ble to the effect in the non-relativistic case. 
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B. Thermalized non-relativistic nuclear particles 

Although massive nuclear particles contribute only a 
negligible portion to the whole energy makeup during 
the period we study, their abundances are crucial in de- 
termining the late time nuclear synthesis. As a result 
of the rapid expansion, the energy density of the radia- 
tion decreases swiftly as in classical treatment, while 
matter density such as that of protons and neutrons de- 
creases much slower as a^^. Additionally, relativistic par- 
ticles and their anti-particles annihilate with each other 
as temperature goes down, and the non-relativistic nu- 
clear particle component, therefore, becomes increasingly 
important. The non-relativistic feature of nuclear par- 
ticles shall preserve as temperature decreases through- 
out the thermalization and into nuclear synthesis. Since 



most of the nuclear particles are fermions, we consider 
only non-relativistic fermion gases without loss of gen- 
erality. We can safely neglect the chemical potential of 
Bose gases since the photon number is not conserved. 
But for non-relativistic Fermi gases, at extremely high 
temperature, such a simplification may no longer be ap- 
propriate. Fortunately, it turns out that the chemical 
potential fails to show up in density and pressure expres- 
sions to the first and second order approximations in our 
result. This not only significantly simplifies the calcula- 
tion, but also justifies the usual negligence of chemical 
potential in the evolution of non-relativistic fermions in 
the early universe 

For massive fermions, the grand canonical partition 
function is 
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We then obtain the energy density in a conventional way is found to be 
U = -^lnE, (23) 



Cl 



P 



U 
V 



(24) 



That is, the density is a function of /3 and volume, while 
/3 and V vary in time. Here the background temperature 
evolution can be well approximated by T ^ t^^/^ dur- 
ing the radiation dominated period. In the time span we 
are studying, /? varies within one order of magnitude. To 
simplify the problem correspondingly, we set the temper- 
ature to be 10"'^^ K and solve the density evolution. The 
mass of general non-relativistic particles in calculation 
is chosen to be the neutron mass, since masses of other 
non-relativistic nuclear particles at that time differ only 
by a coefficient within one order of magnitude. 

The expression for p becomes exceedingly lengthy. Fac- 
ing a huge magnitude range which unavoidably exceeds 
the calculation limit in the computer, to evaluate the ex- 
pression without any numerical reduction is impossible. 
However, after observing the density behavior within a 
small volume range, an asymptotic behavior is found. 
We then analyze the density behavior as volume radius 
increases to four light years. After careful magnitude 
analysis, we find that as long as the volume of the sys- 
tem exceeds 10~^*'T4nit, where Vinit is the initial volume 
of the primordial patch. Such a requirement is thus safely 
guaranteed in the primordial patch before BBN. The 
first-order density behavior relying on the temperature 
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ciT 
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(25) 



where ci equals 1.27 x 10 ^8 in Planck unit, or 3.9 x 
103 J . K5/2/m3. We use the extreme condition of the 
early universe between 0.1 s and 1 s to numerically re- 
duce the expression and find that the result goes back to 
the ideal fermion case, as is shown in Fig.[TJ Applying the 
same method, a similar asymptotic behavior of the first 
order can be found in the pressure- volume relation. We 
notice that strange negative pressure occurs in the trans- 
planckian scale due to the failure of statistical physics. 
However, as the radius of the patch grows to four light 
years, which is the initial value of our calculation, the 
pressure again becomes almost volume- independent. The 
reduced pressure is 



5/2 



5/2 



(26) 



where C2 = 2.60 x lO'' J • /m'^, as is shown in Fig. [2l 

Then the equation of state in the specific time span 
becomes 



wP, 



(27) 



where the state constant w — 1.5 to the first-order ap- 
proximation, as can be seen in Fig. [TJ The recurrence 
back to classical case implies that the gravitational effect 
on non-relativistic gases is quite small and can be safely 
neglected to the first order. We carry on the calculation 



8 




2 4 6 8 10 



FIG. 1: The reduced density /? relation in Planck unit. 




FIG. 2: The reduced pressure /3 relation in Planck unit. 

to the second order and find the density and the pressure 
behave as 

/ 1 \ 5/2 
/ 1 \ 5/2 

where Ci and £2 have numerical values 1.37 * 
10"''° J • K/m^ and 2.06 * 10""^^ J • K/m^ respectively. 
The state constant to the second order evolves in time 

where C3 equals 3.54 * lO^^'"^ K'^/^, with temperature 
T ^ lO^^K. Thus we find the total modification in- 
duced by the self-gravitational effect in the EoS to the 
second order acts as an additional term of the magni- 
tude O(10~*^) and increases over time as c^/T^^^ ^ t^/^. 
Though relatively small, we can see that such a term 



tends to decrease the pressure between non-relativistic 
particles as a manifestation of gravitational attraction. 
Besides, the deviation grows as temperature cools down, 
and finally regains its power in the matter-dominated 
era. The whole analysis serves as a robust test to the 
justification of discarding self-gravity when solving the 
evolution of the non-relativistic nuclear particles during 
the radiation-dominated period. 



IV. REMARKS 

In general relativity, any kind of dynamic evolution, 
if inducing local energy fluctuation, cannot escape the 
gravitational back-reaction, i.e. the change of space time 
geometry itself. From the principle of minimal cou- 
pling, such dramatic feature of general relativity indeed 
does not affect different kinds of interactions significantly 
enough in most cases. Consequently, the theory of gen- 
eral relativity is poorly tested except in a weak-field limit 
within small regions like solar systems and radio pul- 
sars, no bigger than galaxy scale. However, the Fried- 
man equation using Robertson- Walker metric from gen- 
eral relativity is widely utilized in solving the homoge- 
neous and isotropic expansion of the universe. A test 
of general relativity in a large scale of cosmic size then 
proves crucial and necessary in order to justify the foun- 
dation of modern cosmology. One of some promising 
methods is to use particle element abundances [l^l to 
restrict expansion rate before BBN. 

The weak interaction before BBN between neutrons 
and protons should be frequent enough to counter ex- 
pansion in maintaining equilibrium. Accordingly, once 
the expansion speeds up, the weak interaction between 
neutrons and protons will be out of equilibrium at higher 
temperature, causing a larger neutron-to-proton ratio. 
In turn, the abundances of deuterium, helium-3, and 
helium-4 are increased. Knowing that helium-4 is hard to 
destroy during galaxy formation, observations of helium- 
4 abundance can help test the expansion rate and the 
cosmological constant as previous work suggested jTsj . 
Modifications of the Friedman equation in their work can 
be parameterized as 

(^)=-^p[l + x(i^/3pc2)], (31) 

where x corresponds to the modified state constant w in 
our case, see Eq. ((2T|l and Eq. ((30|) . In existing obser- 
vation results, such as those by WMAP 19], Peimbert 
et al. [2(3|, and Steigman [2l|, the accuracy of the con- 
straint of X is of the magnitude 0(0. 1), well beyond the 
self-gravitational effect in our result. Thus, we can ex- 
clude self-gravity as the dominant mechanism behind any 
observable deviation in this method. 

Finally, reexamining the derived formula Eq. (1211) for 
state constant, we compare it with the result in boson 
stars [1], and find that w in both cases increases as the 
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volume of the system expands. The merging similari- 
ties of the self-gravitational effect between two distinct 
circumstances are illuminating to understand the funda- 
mental properties of gravity itself. 

V. CONCLUSION 

In this paper, we try to evaluate the self-gravitational 
effect in Bose and Fermi gases at the early universe and 
set an exact bound on the magnitude of the modifica- 
tion. We study the Bose and Fermi gases before BBN 
synthesis in thermal equilibrium at a temperature in the 
vicinity of 10^^ K, when particles collide relativistically 
with each other frequently enough to maintain thermal- 
ization in spite of the significant cosmic expansion. We 
formulate new grand canonical partition functions based 
on the proof Uhlenbeck and Gropper made on non-ideal 
Bose and Fermi gases. Taking advantage of the severe 
state of the early universe, we reduce the modified den- 
sity and pressure expressions and derive the equations of 
state of relativistic and non-relativistic gases, discovering 
that self-gravity manifests itself in a simple way. In both 
non-relativistic and relativistic cases, we find the devia- 
tion from ideal ones to be small. This justifies the ne- 
glecting of the self-gravitational effect when dealing with 
non-relativistic particles. However, dynamical behavior 
appears when considering the second order effect. Our 
result serves as a robust proof for neglecting self-gravity 
when solving the evolutions of relativistic Bose and Fermi 
gases as well as the thermalized non-relativistic Fermi 
gases. Additionally, we distinguish different effects of 
self-gravity on the relativistic and the non-relativistic 
particles. In thermalized relativistic particles, discrep- 
ancies between boson and fermion statistics become un- 
apparent, and the self- gravitational effect of magnitude 
©(lO"^"*^) decreases with time. This is naturally ex- 
pected as matter becomes more and more dominant in 
the energy content while the energy density of radiation 
undergoes a more rapid dilution. Gravitational attrac- 



tion consequently diminishes much faster in radiation 
components than in matter. While in non-relativistic 
Fermi gases, the self-gravitational effect, though still very 
small of the magnitude O(10~*^), increases over time as 
the temperature cooling down dramatically and finally 
becomes significant in the matter-dominated era when 
both homogeneity and isotropy break down and the EoS 
description is no longer appropriate. As to the evolution 
of the background metric, we find that self-gravity of any 
kind will not strongly affect the conventional result using 
EoS of ideal relativistic and non-relativistic boson and 
fermion statistics. Thus the previously not fully justified 
negligence of self-gravity is theoretically studied from a 
statistical perspective when dealing with cosmic evolu- 
tion before BBN. However, our study is restricted to a 
given temperature range, and the same result may no 
longer be true at the very early stage of the universe as 
can be seen in Eq. (j2ip . At a much earlier epoch, when 
the volume of the universe is small enough with the basic 
assumption of our statistical approach of thermal equilib- 
rium still satisfied, self-gravity may no longer be negligi- 
ble from our result. Thus, further study should be carried 
out at much earlier stages to give a complete evaluation 
on the effects of self-gravity and then the back-reaction. 
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